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Abstract. This paper considers two-dimensional stably stratified steady periodic gravity 
water waves with surface profiles monotonic between crests and troughs. We provide 
sufficient conditions under which such waves are necessarily symmetric. This is done by 
first exploiting some elliptic structure in the governing equations to show that, in certain 
size regimes, a maximum principle holds. This then forms the basis for a method of moving 
planes argument. 



1. Introduction 

One of the characteristic features of ocean waves is their propensity to display stratifi- 
cation — a heterogeneous distribution of density that can resuh from the interplay between 
gravity and the salinity of the water In a typical scenario, one observes a region just below 
the surface where density increases quickly with depth, followed by a larger region extend- 
ing to the ocean bed where it is essentially constant (cf. ifTSl '). Qualitatively, flows with 
heterogeneous density may differ quite strikingly from their homogeneous counterparts. 
For instance, slowly moving stratified traveling waves have a strong tendency to be two- 
dimensional. That is, the motion is identical along any line perpendicular to the direction 
of propagation and the vertical axis 1, 1 8 1 . Remarkably, this columnar behavior will persist 
even with a foreign body, such as an infinitely long cylinder, placed vertically within the 
fluid. 

Unsurprisingly, therefore, stratification has been the subject of a great deal of scholarly 
interest, especially in the oceanography and geophysical fluid dynamics communities. In 
06, 1 7 1 , the author developed an existence theory for two-dimensional stratified steady and 
periodic gravity waves, with or without surface tension. It was shown that, under certain 
unrestrictive assumptions, there is a global continuum of such waves that are classical. This 
was done via a bifurcation argument, with the 1 -parameter family of laminar flows playing 
the role of the trivial solutions. As a consequence of this construction, it was shown that 
each member of the continuum has monotonic surface profile and is symmetric. That is, 
the wave profile possesses a single crest and a single trough per period and is monotonic 
between them. Moreover, both the horizontal velocity and surface are symmetric across 
the crest line, whereas the vertical velocity is antisymmetric. In this work we address the 
complementary question: Under what circumstances are stratified steady water waves with 
monotonic profiles necessarily symmetric? While at first blush it may seem that waves of 
this type are quite special, our main result implies that, in fact, they are the only possible 
solutions within certain size regimes. 
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Let us now briefly review the setup for stratified waves in 116|. Fix a Cartesian coor- 
dinate system so that the x-axis points in the direction of propagation, and the y-axis is 
vertical. We assume that the floor of the ocean is flat and occurs at y = -d. Let y - rjix, t) 
be the free surface. We shall normalize by choosing the axes so that the free surface is 
oscillating around the fine y = 0. As usual we let u - u(x, y, f) and v - v{x, y, t) denote the 
horizontal and vertical velocities, respectively, and let p - p(x, y,t) > Q be the density. 

For water waves, it is appropriate to suppose that the flow is incompressible. Math- 
ematically, this assumption manifests itself as the requirement that the vector field be 
divergence-free for all time, 

(LI) + Vy = 0. 

Taking the fluid to be inviscid, conservation of mass implies that the density of a fluid 
particle remains constant as it follows the flow. This is expressed by the continuity equation 

(L2) p,+UPy, + VPy=Q. 

Conservation of momentum is described by Euler's equations 



(1.3) 



pUt + pUUji + pVUy - —fx, 

PV, + PUV,, + pVVy = -Py - gp. 



where P - P{x, y, f) denotes the pressure and g is the gravitational constant. Here we are 
assuming that gravity is the only external force acting on the fluid. 

On the free surface, the dynamic boundary condition requires the pressure of the fluid 
to match atmospheric pressure, that we shall denote Patm. Thus, 

(1.4) P = Patm, ony^riixj). 

The corresponding condition for the vector field is the so-called kinematic condition 
(L5) v-rjt + urjx, ony-ri{x,t). 

In essence, ( IL5l l states that fluid particles on the free surface remain there as the flow 
develops. Finally, we assume that the ocean bed is impermeable and thus, 

(1.6) v = 0, ony^-d. 

Note that, since we are in the inviscid regime, we do not impose any conditions on u along 
the bed. 

We are considering classical traveling periodic wave solutions {u,v,p,P,T]) to (ILU - 
(|L6b . More precisely, we take this to mean that, for some fixed c > 0, the solution appears 
steady in time and periodic in the jc-direction when observed in a frame that moves with 
constant speed c to the right. The vector field will thus take the form u - u{x - ct, y), 
V = v(x - ct,y), where each of these is L-periodic in the first coordinate. Likewise for 
the scalar quantities: p - p(x - ct,y), P = P(x - ct,y), and // - rj{x - ct), again with 
L-periodicity in the first coordinate. We therefore take moving coordinates 

(x — ct,y) I— > {x,y), 

which eliminates time dependency from the problem. In the moving frame ( |Ll| i-( fT3] l 
become 



(1.7) 



Ujc + Vy - Q 

{U - C)px + VPy - 

p(u - c)Ux + pVUy - -Px 

P(U - C)Vx + PVVy = -Py - gp 
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throughout the fluid domain. Meanwhile, the reformulated boundary conditions are 

{V - (u - c)ri^ ony - ri(x) 
V = on y^-d 

P = P-itm ony^ 7]{x) 

where (u, v,p, P) are taken to be functions of x and 77 is a function of x, and all of them 
are L-periodic in x. 

Waves of this type will be called symmetric provided that there exists xo € R such that 
(u,p,ri) are symmetric over the line x - xq, while v is antisymmetric. We shall call any 
local maximum of 77 a crest, and any local minimum a trough. 

In the event that m = c somewhere in the fluid we say that stagnation has occurred, as in 
the moving frame the fluid appears to be stationary at that point. As will become apparent 
later, stagnation points may produce some degeneracy or instability in the problem. For 
that reason, we shall restrict our attention to the case where u < c throughout the fluid. 

Recall that we have chosen our axes so that rj oscillates around the line y = 0. In other 
words, 

(1.9) f ri(x)dx^O. 

J-L/2 

We shall also denote 

y/max := max //(x) + d, //min := min r](x) + d. 

These are the maximum and minimum distances between the surface and the bed, respec- 
tively. 

Observe that, by conservation of mass and incompressibility, p is transported and the 
vector field is divergence free. Therefore we may introduce a (relative) pseudo-stream 
function t// = t//(x,y), defined uniquely up to a constant by: 

i/^x = - ^/pv, 4iy = ^lp{u - c). 

Here we have the addition of a p term to the typical definition of the stream function for 
an incompressible fluid. This neatly captures the inertial effects of the heterogeneity of the 
flow ifTSl . The particular choice of yfp is merely an algebraic nicety. 

It is a straightforward calculation to check that ij/ is indeed a (relative) stream function in 
the usual sense, i.e. its gradient is orthogonal to the vector field in the moving frame at each 
point in the fluid domain. As usual, we shall refer to the level sets of ^ as the streamlines 
of the flow. Observe that, in assuming m < c, we have guarantee that the streamlines are not 
closed. For definiteness we choose if/ = Q on the free boundary, so that 41 = p^ony - -d, 
where po is the (relative) pseudo-volumetric mass flux, 

nnix) 

(1.10) Po := ylp(x,y)[u(x,y) - c]dy. 

J-d 

It is easy to check that po is well-defined, i.e., the integral on the right-hand side above is 
independent of x. Examining dl.lOl ), it is clear that po describes the amount of fluid flowing 
through any vertical line extending from the bed to the free surface and with respect to the 
transformed vector field ( -\/p(u - c), yfpv). 

Since p is transported, it must be constant on the streamlines. Abusing notation we may 
therefore let p e C^({p^), 0]; R^) be given such that 



(1.11) 



p{x,y) ^ p(-tf/(x,y)) 
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throughout the fluid. When there is risk of confusion, we shall refer to the p occurring on 
the right-hand side above as the streamline density function. We shall focus our attention 
on the case where the density is nondecreasing as depth increases. This assumption is 
physically well-motived, since it is a prerequisite for hydrodynamic stability. Note that the 
level set -if/ - po corresponds to the flat bed, and the set where -i// - corresponds to 
the free surface. For that reason, we shall say the flow is stably stratified provided that the 
streamline density function is nonincreasing. 
By Bernoulli's theorem, the quantity 

P 
2 

is constant along streamlines. Then, under the assumption that u < c throughout the fluid, 
there exists a function /? e C'([0, |;?ol]; K-) such that 

dE 



. c)2 + v2 



) + 8py, 



(1.12) 



ix,y) = -j3(il,(x,y)). 



For want of a better name we shall refer to /3 as the Bernoulli function corresponding to the 
flow. Physically it describes the variation of specific energy as a function of the streamlines. 
It is worth noting that when p is a constant, /3 reduces to the vorticity function. 
We now define several quantities describing various properties of the wave. Let 

i>2 



:= eM? : -d <y < t](x)] 



denote the fluid domain (in the moving frame) and, for each € N, domain Q c 
introduce the space 



we 



Put 



ao :-- 



{/ e C*(Q) : fix + L,y) ^ f(x,y), ^x,y G q} . 



rnin ^fp{c - u), 



Aq :- rnax \fp{c - u). 



The first of these gauges how far away from stagnation the flow is, while the second gives 
a bound on m , the motion of the fluid counter to the direction of propagation. Next, we let 



(1.13) M ;= max< 





V 




max 




, max 






u — c 





\ u — C / U - I 



rnax 



1 



■\/p(u - c) ' u — c 



M describes (in a sense that will be made clear later) the variation of dy/dx following along 
a streamline. In particular, M = only for the family of laminar flows, these being shear 
flows where the streamlines are parallel to the bed and 77 = 0. The motivation here will 
become significantly more transparent in the next section, in particular, see (12.7b and the 
accompanying discussion. 

Our main result is the following. 

Main Theorem. Consider a stably stratified steady periodic wave train propagating at 
fixed speed c over a flat bed at y — —d with relative pseudo-mass flux po < 0. Let 
p e C^([po, 0]; R"*") and fi e C'([0, \pq\\) be the streamline density function and Bernoulli 
function associated with the flow, and let (u, v) e Cp^^(D^) x Cp^X^,;) be the vector fleld. 
Assume that the wave proflle y — ri(x) is monotonic between crests and troughs with period 
L and that u — c < 0. Each of the following is a sufficient condition for the wave to be 
symmetric. 

2 



(SI) 



'7max ^ 



sup/3' + gri^^ sup (p'Y < n 



0< i< Ipol; 
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or 

(S2) supyS' + CTmax sup (p")^ < cxp (- min{L, ;?min)); 



or 
(S3) 
where 
and 



^0^1 a/S + ^2 ^ + g sup Ip'l < exp (-fl„ '^^ min{L, |pol}), 
ei < ^, ^2 < ao. 

ei := 3CTmaxSUp|p'| , 



62 := AMAl + IM^Al + 3Ao (sup |j8| + ^max sup |p'|) 

+ IM^a^^Al + M^flQiAf) + MflQ^Af) (sup + ^?7„,x sup |p'|) . 



Remark. The first condition ( ISII ) is a direct generalization of the main resuh of |4|, while 
the second follows from a similar argument but takes a slightly different view. We mention 
that, in particular, the condition /?', p" < immediately implies both of conditions dSlb - 
( IS2b . Note that this is a physically natural assumption on p for ocean waves fr2|. For/?' or 
p" small but positive, the Main Theorem implies that the wave is symmetric provided that 
the maximum elevation of the profile and/or the period are sufficiently small. 

On the other hand, condition ( IS3I ) gives a criteria based on the distance from stagnation, 
the degree to which the flow is laminar, as well the magnitudes of the volumetric mass 
flux, Bernoulli function, density variation and the period. Generally, it states that solu- 
tions exhibiting small fi and p' and which are sufficiently far from stagnant are necessarily 
symmetric. 

Let us now discuss the history of this problem. The basic machinery we employ is the 
moving planes method, which was originated by Alexandroff 1 1 1 . Since then, this argument 
has been used extensively by a number of authors, notably Gidas, Nirenberg and Ni 1 8 1 and 
Serrin |14|. The first significant application of a moving planes argument to the study of 
water waves was by Craig and Sternberg |5|, who considered the case of solitary waves in 
the irrotational setting. 

The direct inspiration for our work is a number of recent breakthroughs in the study of 
symmetry for two-dimensional constant density rotational gravity water waves, beginning 
with a result of Constantin and Escher [4] analogous to condition ( ISll l of our main theorem. 
This was a generalization of an earlier result of Okamoto and Shoji, who showed in ifTTI 
that symmetry must occur when the profile has a single crest and is monotonically decreas- 
ing. Soon thereafter, Hur, building on the work of Garabedian |7| andToland |15| for irro- 
tational waves, proved symmetry for arbitrary vorticity functions under the assumption that 
every streamline attains a minimum on the trough line |9|. While this last condition was 
not entirely physical, she also introduced the idea of working within the Dubreil-Jacotin 
formulation, which proved to be extremely convenient. As we shall see, the main advan- 
tage of Dubreil-Jacotin is that it demotes the vorticity function from a semilinear forcing 
term to the coefficient of a first-order term in a quasilinear equation, thereby nullifying its 
importance for symmetry. Refining Hur's argument, Constantin, Ehrnstrom and Wahlen 
were able to significantly improve upon the result in 14]. In (|3J, they showed that every 
two-dimensional rotational steady gravity wave with a monotonic profile is symmetric. In 
other words, symmetry follows without any assumptions on the vorticity function. 
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Work on the heterogeneous case has been comparatively sparse. The most complete 
results are due to Maia ifTOl . who considered solitary internal waves confined to an infinite 
horizontal strip. It was assumed that the waves were quiescent at +00 and connected to 
a laminar flow at -00, where a piece-wise continuous density distribution was prescribed. 
Maia proved that, in this setting, waves of elevation are symmetric. Her argument, however, 
relies heavily on both the geometry of the domain and the type of boundary data at ±00, 
and therefore does not extend to the free surface and periodic case which we consider. 

The structure of our argument is as follows. In ^2.1l we begin by providing two scalar 
reformulations of ( ll.7b - (ll.8b : these are the natural generalizations to heterogeneous flows 
of the equations studied in [3J and [4|. In ^2.21 we then derive estimates on the given 
quantities that ensure solutions of these equations satisfy a maximum principle. This sets 
the stage for the moving planes argument that we provide in ^ 

2. Preliminaries 

2.1. Two reformulations. It will be convenient to replace the system (ll.7b -( ll.8l l with an 
equivalent scalar PDE for the stream function. Differentiating the expression for E and 
appealing to the definition of if/, one can prove the following identity, known as Yih's 
equation or the Yih-Long equation ITSl . 

dE ^ dp 

dip dip 

Rewriting the above expression in terms of the Bernoulli function jS and streamUne density 
function yields the following 

(2.1) -l3{ip)^Ailf-gyp'(-ilf). 

Moreover, evaluating Bernoulli's theorem on the free surface s 0, we find 

(2.2) IViAl^ + 2gp{-iP) (rjix) + d)^Q, ony ^ rj{x) 

where the constant Q :- 2{E\,j - P^tm + gd). Note that Q gives roughly the energy density 
along the free surface of the fluid. Together with the fact that \p - -po on y - -d and 
(/r = on the free surface, ( |2.1b -( |272] l provide a complete reformulation of the problem as 
a semilinear elliptic equation. 

It will also prove useful to consider the semi-Lagrangian (or Dubreil-Jacotin) formu- 
lation, wherein we trade some additional nonlinearity in exchange for a fixed the domain. 
Consider the alternate coordinate system, (q, p), where 

q^x, p^-tp(x,y). 

Then, under the transformation {x,y) i-> {q,p), the closed fluid domain is mapped to the 
rectangle 

R := {(q,p) e : -L/2 <q< L/2, po < p < 0). 
The purpose of the minus sign is simply to flip the rectangle so that the free surface will 
correspond to the top of R, while the flat bed will mapped to the bottom. Given this, we 
shall denote 

T := {iq,p) eR:p=0}, B := {iq,p) e R : p ^ po}. 
Note that, in light of (II. lib and (11.12b . we have that 

P^Pi-p), p^p(p). 
Moreover, the assumption that the streamUne density function is nonincreasing becomes 

(2.3) pp < 0. 
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Next, following the ideas of Dubreil-Jacotin, define 

(2.4) h(q,p):^y + d 

which gives the height above the flat bottom on the streamline corresponding to p and at 
X - q. We calculate; 

1 h„ 

(2.5) «Av = -^, -A. = 7^, 

hp hp 

and 

(2.6) dq = + hqdy, dp = hpdy. 

This is the motivation for the definition of M in (11.131 1. Using (12.5b to write (12.61 1 in terms 
of the Eulerian quantities (u, v), we see that 

(2.7) M = ||/;^|lc,(^). 

For a laminar flow, of course, we have hq, hqq, hqp = 0. It is in this sense that M describes 
the distance of the flow from laminar 

Observe that ( I2.5l l also implies hp > 0, because we have stipulated that u < c throughout 
the fluid. Yih's equation ( 12. Il l becomes the following 

(2.8) - hl/3{-p) = ( 1 + hl)hpp + hqqh] - Ihqhphpq - g(h - d)hlp p 

where we have used ( 12.4b to write y - h - d. Recall, however, that we have normalized rj 
so that it has mean zero. Taking the mean of ( 12.4b along T, we obtain 

nL/2 

(2.9) d = d(h) = f h{q,Q)dq. 

J -L/2 

That is, the average depth d must be viewed as a linear operator acting on h. Namely, it is 
the average value of h over T. Where there is no risk of confusion, we shall suppress this 
dependency and simply write d. 

Next consider the boundaries of the transformed domain. On the bed we must have, by 
the definition of h, that 

(2.10) h = 0, onB. 

In the new coordinates, the definition of Q given in ( 12.21 ) becomes the requirement 

(2.11) l+hl + hl (2gph - 6) = 0, on T. 

Altogether, then, combining ( 12.8b . ( 12.101 ) and ( 12.111 ) we have that the fully reformulated 
problem is the following. Find {h, Q) € C^^^{R) x R satisfying 

[■ (1 + hl)hpp + hqqh], - Ihqhphpq - g(h - d(h))hlpp = -hl/3(-p) pQ< p <0 

(2.12) I l+hl + hl(2gph-Q)^0 p^O, 

[ h=0 Po, 



where hp > 0. Here p e C^([po,0];R^) and /3 e C'([0, |/?ol]) are given function with 
Pp<0. 
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2.2. Maximum principle. In anticipation of our argument in the next section, we now 
attempt to find an elliptic problem satisfied by the difference of two solutions of the height 
equation. Fix h, h solutions to ( I2.12l i with the same value of Q and suppose that d{h) - diji). 
Define 

X := {\ + h^^hfd], + h-'dl-2hqhfdpdq 

+ {hqqihp + hp) - Ihqhpq + fi(-p){h], + h,}ip + ypp^hfdp 

(2.13) + i^ppihq + hq) - Ihphpq) hfO, 

-gpp(h - d(h))(hl + hphp + hl)hfdp 
-gPp- 

Note that the first three lines of ( 12.131 ) are the same as the analogous operator found in lO, 
except that we have divided through by h^p\ the stratification manifets itself only in the last 
two lines. Also, since pp < 0, the sign of the zeroth-order term "goes the wrong way" 
in the sense that, even taking for granted that X. is elliptic, the maximum principle does 
not generally apply. The goal of this section is to find suitable restrictions on h so that a 
maximum principle like result will hold. Our main tool for this is a number of results due 
to Berestycki, Nirenberg and Varadhan We begin by summarizing the most relevant of 
these. 

Definition 1. Let X. be an elliptic second-order linear diff'erential operator with bounded 
continuous coefficients. The principal eigenvalue Ai - Ai{X.) is defined as follows 

Ai sup{^ € R : 30 > in i? satisfying (£ + A)4> < 0). 

The fact that Ai is well-defined, i.e., that the supremum above exists, is a classical result. 
Our interest in Ai stems from the following theorem. 

Tlieorem 1. (BNV, JJl) Let X. be an elliptic second-order linear differential operator with 
bounded continuous coefficients. Then X has a ( refined) maximal principle if and only if 
AiU) > 0. 

The meaning of "refined" above is made precise in [2 |, but for our purposes it can be 
safely said to coincide with the traditional definition of the (weak) maximum principle. 

In light of Theorem [U our objective is clear: We seek to prove lower bounds on the 
principal eigenvalue of the operator in ( 12.13b . The most basic such estimate is given by the 
following result of Protter and Weinberger lfT3l . 

Tlieorem 2. ( PW, II13I ) Let XL be a general second order linear differential operator with 
bounded coefficients. Then, if<p e C^{R) n C'(/?) is a strictly positive function on R, 

Using this with together with a construction in (]2l we can prove a very rough but explicit 
bound on Ai in the case where there is no zeroth-order term. 

Corollary 1. Let M :— ajjdidj + bjdi be a second-order elliptic linear differential operator 
with bounded continuous coefficients. Let ao denote the lower ellipticity constant and put 
b^ :— Yj b]- Let cr > solve the equation 

aocr^ - bcr - b - I. 
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Then 

(2.14) Ai(M) > exp(-o-min{L, \pq\)). 

Proof. In 121 it was shown that there exists a function mq satisfying Muq = - 1 in m > 
in R and u\aR - (potentially in some weak sense.) If we further examines the actual 
construction of mq, we find that < mq < e°"™"(^'l'^oll, with cr as above (see §3 of [2J, 
particularly p. 61). Then (12.14b follows directly from Theorem|2]with (p = uq. □ 

Finally, for reference we repackage two important estimates from |]2l regarding the 
dependence of Ai on the magnitudes of the first- and zeroth-order terms. 

Proposition 1. ( BNV) Let M := aijdidj + bidi and M' :- aijdidj + b'.d; be a second- 
order elliptic differential operators with smooth coefficients. Let ao be the lower ellipticity 
constant and put b^ ;= 2i Suppose that 6 is chosen so that 

2 - bf <6^< bao, 

i 

then 

(2.15) M(M') > Ai{M) - J— 6. 

V ao 

Proof. This is essentially estimate (5.2) of f2] (used to prove Proposition 5. 1 of that paper) 
with c = 0. □ 

Proposition 2. ( BNV) In its dependence on the zeroth-order term, Ai is Lipschitz contin- 
uous in the L°°-norm with Lipschitz constant L 

That said, we are now prepared to prove the key lemma of this section. 

Lemma 1. (Properties of £) Let h, h e Cll^{R), Qe^be given such that d(h) - d(h) 
and (h, Q), (h, Q) are solutions of ( |2.12| l with sup^/;^' - sup^/i^', inf^/z^' - inf^/i^' . 
Then the operator XL defined in ( 12.131) satisfies the following. 

(a) Uh-h)^Q. 

(b) £, is a uniformly elliptic differential operator with lower ellipticity constant ao — 
mfjhj,\ 

(c) Let /li(-£) denote the principle eigenvalue of X. as in ([T]). Define 

b := sup!.[h-/h^g(hp + hp) - Ihfhqhpq +l}{-p)hf{hl + hphp + hj) 

(2.16) - hp'gppih - d{h))(hl + hphp + hl)f 

+ [hp^hppihg + hg) - 2hp%gf 

and let cr be as in Corollary\l\ Then A\{Ji) is strictly positive provided that 

(2.17) g-'"™"'^' > sup Igp^l. 

Alternatively, if we recall Aq — sup^/z^' and define the constants 

(2.18) 5i := 3|/?ol sup \gpp\, 

[po,0] 



n2 1 2 
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and 

62 := sup I [h-/hqq{hp + hp) - Ikflijipq +l3{-p)hf{hl + hpJip + li])^ 
(2.19) " , 

+ ^hp^Jippihg + Jig) - Ihp^hpqj \ , 



then is strictly positive provided that 



(2.20) ^1 < ^2 < ao, 



and 

(2.21) Aq5i J— + 62 J— + sup \gpp\ < exp (-Oq min{L, \pq\]). 

y ao y ao [p^^] ' 

Proof. A tedious but easy calculation readily confirms that £Xh - Ii) - 0, proving (a). 
Similarly, the second-order terms of X. are identical to those of the differential operator in 
the height equation which has been proven to be uniformly elliptic lfT6ll . Calculating the 
precise lower ellipticity constant is a simple exercise that we omit for brevity. 

The interesting part is, of course, proving (c), which is done as follows. By Proposition 
|2]we know that Ai is Lipschitz in its dependence on the zeroth-order coefficients with 
Lipschitz constant 1 . That is, if we denote as Xo the differential operator found by dropping 
the zeroth-order term in then 

(2.22) AM) > MUo)- sup \gpp\. 

In order to show Ai{X.) is positive, therefore, we must produce lower bounds for Xo. This 
can be achieved by using Corollary [T] directly, proving Ai{X.) is positive under condition 
dTTTl i. 

To arrive at the more involved conditions ( I2.20I )-( |Z2TI ) we instead begin by chopping 
off the first-order terms involving pp using Proposition[T] Let denote the operator found 
by setting pp = Oin X.- Let b be as in ( 12.161 ) and put 

5 sup \gpph-/ih - d{h)){hl + hphp + hl)\ . 

R 

Under the assumption that < bao we have by (I2.15l i: 

(2.23) Ad-Co) > AM - J -S. 

V ao 

So that we can better understand the above statement we now endeavor to control 6. We 
can bound h - d(h) hy ?] + d, where 77 is the free surface corresponding to h. Since we have 
stipulated that supj hp = sup^/j^, 

6 < 3Ao sup \gpp\ sup \fi + d\. 

Ipo.O] T 

AsS < b, the assumption < bao can be satisfied by requiring that 6 < ao, or equivalently 

3 sup \gpp\ sup \ri + d\< flo/Ao. 

[po,0] T 

We remark that because ao < Ao, the left-hand side above is necessarily less than one. It 
warrants notice that this implies there is a definite bound on either the magnitude of the 
density variation, or alternatively on the amplitude of the wave. We can eliminate fj + d if 
we are wilHng to further estimate ij + d(h) < \po\ sup^/i^ = \po\a^^^ . Plugging this into the 
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above inequality reveals that it is sufficient to take 6i = 3\po\ sup \gpp\ < a^/Ao in order to 
ensure estimate (12.231 1 holds. 

Using this information in concert with ( 12.221 ) we obtain 

[T 

AiU) > Aii£i) - sup \gpp\ - J—S 

V ao 

V flo 

Next, let us denote by X.2 the operator consisting of the second-order terms of X.- Ap- 
plying Corollary[T]to X2 yields 

The final step is to estimate the difference between /li(-Ci) and Ai(X.2)- By the same argu- 
ment as (12.23b we get 

(2.24) MUi) > AM - J^Si, 

V ao 

since we have 62 < ao by ( 12.201 ). 

Incorporating the last estimate with ( |2.22| l-( |Z23] l. yields the condition in ( 12.211 ). Our 
two assumptions on the relative size of aq and (61, 62) — made when deriving ( 12.231 ) and 
( 12.24b — are collected in ( 12.201 ). This completes the proof of (c) and the lemma. □ 

3. Symmetry 

We are now prepared to prove the main theorem. This shall be done by employing 
an adapted moving plane method (cf. Il3ll4li8 , 14,1.) The necessary maximum principle 
properties will follow from the results of the previous section. 

Proof of Main Theorem. First consider the case when dSlI ) holds. We shall structure our 
argument on those of f4l. Without loss of generality, suppose that the trough of the wave 
train occurs at x = ±L/2. Also, since we will use the Euler formulation exclusively for 
( ISlb -(IS2b. we shall make a change of variables so that the bed occurs at y = and the free 
surface is y = rfix) + d. This has the convenient effect of making y non-negative. 
By assumption, we may choose 6 e (0, 1/2) satisfying 

(3.1) l3'{s)T^l.,, + gp"(-s)Til.,,<nHl-2Sf , < i < |pol, 

We may therefore define a function a e C^([0, /ymaxl) by 



a{y) :- -25)-^ + dV 0<y<rimax- 

One can readily verify that a has the following two properties 

ayy _ n\\-25f 

~ 2 ~ 2 ' 

^ '/max '/max 



(3.2) < sin(;r(5) < a < 1, 



For each A e (-L/2, 0] consider the truncated fluid domain 

D,x := £> n {(x,y) : -L/2 < x < A}. 

Since the surface profile is monotonically decreasing with the trough occurring at x = 
-L/2, for A sufficiently near -L/2, the transformation (x, y) {2A - x, y) maps into D. 
For such A we define the reflected domain 

Da U {(2^ - x,y) : (x,y) e D,} c D. 
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Incrementing A, we eventually reach a point where is no longer contained in D. There- 
fore, if we put 

Ao := inf [a e (-L/2, 0] : c d} , 

We consider three possibilities, the first being that Aq = 0. If this fails to be true, then 
Ao < and we further consider the cases where Aq occurs at a crest, and where Aq occurs 
strictly to the left of a crest. 

Case 1: Ao = 0. For each {x,y) e D*^ -: D., put 

ifr{-x,y)- tfr(x,y) 

w(x,y) — , 

a(y) 

where t// is the pseudo-stream function. Proving symmetry of the wave is equivalent to 
showing that w = in f2. Note that t//(-x,y) > il/{x,y) on the free surface, by construction, 
and i^(-jic, y) - ifr{x, y) on the bed and under the trough line, by periodicity. In other words, 
w > on 5Q. This suggests that, in order to prove w vanishes identically, we need to show 
that w satisfies a maximum principle. Naturally, this will have to come from ( 12.11 ). 
With that in mind, we compute 

A o'^>';3 A(aw) a„ . 
Aw + 2 — OyW = —w, m n. 

a a a 

Since a(y)w(x, y) = i//(-x, y) - ^{x, y), from ( 12. It and (13. 2t we conclude 

or,. 

Aw + 2—dyW + cw < 0, in Q., 

a ' 

where 

- _ ;r2 
c -.^p-gyp' - ^— , 

'7max 

and, for all (x, y) e Q, 

P(x,y) ■.^/3(4r(-x,y))-p(if,{x,y)), '^(x,y) p'i-t/,{-x,y)) - p'{-if,(x,y)). 

But, by the mean value theorem, for each ix,y) e Q., there exist si = si{x,y) € (0, |pol), 
S2 = S2(x,y) e (/?o,0), such that'^(x,y) = /3'(si), p'(x,y) = -p"(s2)- Since < y < Ty^ax, 
we have 

c(x,y) < j3'(si) + gi]^,, (p"(s2)y - — , '^(x,y) e Q 

which is non-positive by condition dSlI ). We are therefore able to apply a minimum princi- 
ple and conclude either w > 0, or w = in O. 

Let us denote the corner point C ;= i-L/2, ri{-L/2)). Then, by periodicity of the wave 
profile, we must have 

W(C) = W,.(C) = Wyy(C) = W„(C) = 0. 

Likewise, diff^erentiating the relation = ony = ri(x), we find 

i/r^ -I- 1/',./?' =0, on y = 77(.x). 

Since the trough occurs at x = -L/2, we have ri'(-L/2) = iAx(0 = 0. It follows that 
w,(C) = 0. 

Finally, consider w^y(C). From the nonlinear boundary condition ( 12.2b we compute 

+ ^xyJ]') + (ffyi^xy + ^yyT]') - gp' i-l//)il/x(r](x) + d) + gp{-^)r]' = 0, OU J = r]{x). 
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Evaluating this at C, we conclude if/xy{C)4'yiC) - 0. But, since ij/v < 0, this can only be the 
case provided if/xx vanishes at C. In light of this, and the fact that Wx{C) - 0, it follows that 

Wtv - -2 ■ -Wv = 0, at C. 

a a 

Thus all derivatives of w of order up to two vanish at C. This violates Serrin's edge lemma 
(cf. ||6|), unless w = in O, as desired. 

Case 2: The crest occurs at some point {Aq, rjiAo)), with —L/2 < Aq < 0, and the line 
X = Aq is normal to the free surface there. Put Ai :— 2Aq + L/2, A2 '.— Aq + L/2. We wish 
to consider the domain that results from reflecting D^, over the line x - Ai, and the set 
D \ D 1, over the line x - /I2. More precisely, let us define 



and 



ri(2A() - x) Aq < X < Ai, 
7/(2/12 - x) Ai < X < A2, 



Q := {(x,y) e : /lo < Jc < /I2, -d <y < fjix)}. 

As the crest occurs at (Aq, t](A())), the profile must be non-increasing on the interval x e 
[Ao, L/2]. We therefore observe that Q c D. Let us further define, for all (x,y) e Q, 

i//(x,y) - ip(2Ao - x,y) 



w(x,y) :- 



-W) ^o<x<Au 

Hx,y)-^(2A2-x,y) ^ ^ _^ < ^ 
aiy) 



Suppose that we can prove w s in Q. Explicitly, this means ifrix, y) - if/(2Ao - x, y) in 
Q,, for Aq < X < Ai, and i//(x,y) = 1^(2/12 - x,y) in Q. for Ai < x < A2. Using the fact 
that ifry < 0, we can diff'erentiate the relation (//(x, ri(x)) = to conclude t](x) = fjix), for 
X e (Ao, A2). But this is impossible, as the free surface is assumed to be monotonic. Thus 
w s in Q implies Aq = 0, and so we are in the first case. 

In order to finish the proof of Case 2, we now argue that w vanishes identically. Let C 
denote the crest point (Ao, ?7(/lo)) which lies at the upper left-hand corner of Q. Arguing 
as before, it is easy to show that w > on the upper boundary and w = on the side and 
lower boundary of Q. Likewise, simply computing as before, we show that under assump- 
tion dSlI ). w is a supersolution in Q with respect to the operator A + 2a^^aydy + c, where 
c < 0. Similarly, we check that all derivatives up to order two of w vanish at C. Then, since 
by definition w(C) - 0, Serrin's edge lemma produces a contradiction unless w s in Q. 
Together with the argument of the previous paragraph, this completes the second case. 

Case 3: is internally tangent to the upper boundary of D at some point (^1,77(^1)). 
Clearly, if A^ - 0, then the argument of Case 1 applies and there is nothing additional 
to prove. Without loss of generality, then, suppose Aq < 0. Also, we note that ^ must 
necessarily lie in some interval where the free surface is decreasing. This is a simple 
consequence of the fact that, left of the line x = Aq, the boundary of D^^ is negatively 
sloped. We are therefore justified in defining Ai, A2, 1), O and w as in Case 2. 

We know that i/^ = on the free surface and i/r > in the interior of the fluid. It follows 
that w > on the upper boundary of D.. As before, periodicity ensures that \jj vanishes on 
the lateral boundaries of Q, by construction. Finally, we know that i/r = po any - -d, 
hence w vanishes identically on the lower boundary of Q.. 
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As in the previous cases, a simple calculation shows that w is a supersolution in Q with 
respect to a linear elUptic operator of the form A + c 'o-yC?,. + c, where c < when dSlb 
holds. It follows that w > in Q, or w = in Q. Recall that we have akeady seen that the 
latter of these two possibilities coincides with Case 1 and therefore implies symmetry. 

Seeking a contradiction, suppose instead that w > in Q. Let C := (^1,77(^1)) denote 
the contact point and note that it is in the interior of the upper boundary of Q.. We claim 
further that w(C) = 0. To see why this is the case, observe that by definition of ^, we have 
Tji^) - fii^)- Therefore, in light of the fact that ^ - Qon the free surface, 

iPi^, jj(2Ao - ^)) = = = if,i2Ao - ^, Ji(2Ao - ^)). 

Thus w attains a minimum at C. Since Df^ is internally tangent to D at C, we must have 
that the interior sphere condition is satisfied there. The hypotheses for Hopf's boundary 
lemma (cf. |i6J) now completely verified, we infer that 

dw 

where v is the normal vector at C. 

We seek to produce a contradiction by showing that, in fact, the normal derivative must 
vanish. Towards that end let us consider ^0 2A{) - that is, the reflection of ^\ over 
the line x - Aq. By definition of ^1, we must have that 77(^1) = 77(^0), and 77' (^1) = -ri'i^o). 
Differentiating the relation i^{x, ri(x)) - and taking these observations into account, we 
arrive at the following identity 

(fo.'Kft)) "Ay C 

Here we again made use of the fact that i/r,. = m - c < 0. Also, since p is constant on the 
free surface, the nonlinear boundary condition ( 12.2b implies 



< 0, 

c 



(3.3) - 



(|V^|2 + ^p(-<A)(/7 + d)) = e = (|ViA|2 + gp(-,p)(Ti + d)) 

(ft 



whence 



(3.4) |V.A(^o, 77(^0))!' = |V.A(C)|l 

From (|3.3b -( U4l i we conclude further that 

^Mh lim = -«A.v(C), iAv(^o, '7(^0)) = 'Aj(C). 

Let the normal vector at C be v = (vi, V2). Then, 
dw 



dv 



iAx(C) + <A.v(^o, 77(^0)) iAv(C)-iAv(^o, 77(^0)) «v(77(^i)) 

= Vi h Vo— '■ Vo— -W{C), 

c aijji^i)) a(77(^,)) ff(77(^i))2 



which is zero, by the above considerations. We conclude that w = in Q. This completes 
the proof that condition ( IS lb implies symmetry. 

To prove ( |S2b is sufficient, we note that the only use of dSlI ) in the preceding argument 
was to ensure that a minimum principle holds for w. If we take a = 1, then w will be a 
solution of a linear elliptic operator of the form A + c. Since the ellipticity constant here is 
one, Lemma[T]and Proposition|2]together show that the requirement 

(3.5) supc+ < exp(-min{L, 77inin)) 

n 

implies w has a minimum principle, and hence the wave is symmetric. Condition ( |S2b is 
chosen precisely so that ( 13. 5b holds. 



SYMMETRY OF STRATIFIED WATER WAVES 



15 



The last condition, ( IS3b . will require more effort. In particular, we shall work with the 
height equation formulation (12.12b and make use of the operator X. In doing this we are 
closely following the path set forth by (3\. 

Let h be the height function for a steady periodic gravity water wave as in the statement. 
If h is flat there is nothing to proves, so without loss of generality we suppose that a crest 
occurs in the interval [0, L/2) and a trough at ^ = -L/2. This is permissible, as (12.12b is 
invariant under the transformation q i-> -q. 

For each A € (-L/2, 0), {q, p) e [A, 2A + L/2] x [po, 0], put 

q'^ :— 2A- q, w(q, p; A) :— h(q, p) — h{q^, p). 
Then, by construction, for A e [-L/2, L/2], (q, p) e [A, 2A + L/2] x [0, po], 

(3.6) w(A,p;A)^h(A,p)-h(q\p) = 0, 
and 

(3.7) w(q, po; A) = h(q, po) - hiq-^, po) = 0. 
Also we observe that 

(3.8) d(w(q,p;A)) ^ j- [h(q,0) - h(q\0))dq = 0, 

owing to the periodicity of h. 

Fix q. Then, by monotonicity of the profile, w(q, Q;A)>Q for A near -L/2. Put 

^0 sup [A e [-L/2, 0] : w(q, 0; A) > 0, for q e (A, 2A + L/2)] . 

There are two cases to consider 

Case 4: Aq = 0. Denote Q := (0, L/2) x ipo,0) and, for concision, fix A = Aq for w. 
Identity ( 13.7b implies that w = on the set dQb ■- {(q,Po) '■ < q < L/2]. Likewise, 
by periodicity of h we have w{L/2,p) - h{L/2,p) - h{-L/2,p) - 0, for all po < p < 0. 
Thus w vanishes identically on the right boundary, dQ.,- ■- {{L/2,p) : pQ < p < Q]. Taking 
A and replacing it with Aq = in (13.6b . we obtain the same result on the left boundary, 
dQ.1 :- {(0,p) : Po < p < 0]. Finally, on the top, dQ., :- {(q,0) : < q < L/2], we have 
w > by the definition of Ao. Combining then, we have w > on c?0. 

Suppose for now that w has a determined sign in some region O' c Q whose boundary 
includes the corner point (L/2, 0) and which is blunt in the sense that the Serrin edge 
lemma holds there (see, for instance. Definition E3 of |6|). By (13.6b and periodicity, 
w(L/2, 0; Aq) - Wp(L/2, 0; Ao) = Wq(L/2, 0; Ao) - 0. Using these facts we can differen- 
tiate ( 12.1 lb in q and evaluate at the corner to 

2hphp,{2gph-Q) = 0, at (L/2, 0), 

which implies hpq — and thus Wp^ — vanishes at (L/2, 0). 

As w is a solution for X in Q.' with determined sign, applying Serrin's edge lemma and 
the arguments of the previous paragraph we see that at least one of w^^ and Wpp must be 
nonvanishing at (L/2, 0). This leads to a contradiction, as w(q, p; Aq) = h{q, p) - h{-q, p), 
and thus 

w,,(L/2,0) = V(L/2,0)-/i,,(-L/2,0) 1 
Wpp(L/2,0) = hpp(L/2,0) - hpp(-L/2,0) j 
We conclude that no such region Q' can exist. 

Case 5: Aq e {-L/2, 0). In this case, there also exists a first point qo e (Ao, 2Ao + L/2] for 
which ^(^0, 0; Ao) - 0, and w{q, 0; Ao) > for Ao < q < qo This will be attained at the 
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point where the graphs of the functions q i-> h{q, 0), q i-> h{q 0) are tangent to each other. 

First observe that, by definition of /Iq, we cannot have w{q, 0; Aq) > for q > Aq. We 
may continuously extend w by letting 



L L 

2Aq + -, Aa + - 



x[po,0]. 



w{q, p; Aq) :- h(q, p) - h{q ° + L, p), for (q, p) e 
Redefine 

Q:=^.lo, Ao + ^x{po,Q), 

and, as in f3l, periodicity implies w € C^(0). Also, if 2A + Ljl lies to the left of the crest 
line, w{q, 0; /i) > holds for A < q < 2A + Ljl. Thus 2Aq + L/2 must be to the right 
of (or coinciding with) the crest line. Monotonicity of the profile then ensures that h is 
decreasing for q e (2A() + L/2, L/2). This implies that on the top 

LI 



w{q, 0; Aq) > 0, e 

On the boundary, we therefore have the following sign information 

i w{Ao,p;Ao) = w(/lo + f = 0, for /? e [/?o, 0], 

1 w{q,pQ;AQ) - and 0; /Iq) > 0, for ^ e /to + | j . 

Suppose there exists a region Q" c Q. where w > in Q", (^q, 0) e 5Q" and 50" is 
smooth at (g'o,0). Since the graphs of ^ i-» h(q,Q), q i-> h{q^°,Q) are tangent at the point 
(^■0,0) this point, Wq{q{),0; Aq) - 0. In particular, hq(qo,Q) - -hg(q'*'',Q). Using this to 
evaluate jZTT] ) at (^o, 0) and (q^^Q) we find hp{qo, 0) = hpiq^", 0). Thus Vw(qQ, 0; Aq) = 0. 
However this contradicts the Hopf lemma (which is applicable, despite the adverse sign of 
the zeroth-order term of X, because of the predetermined sign of w in Q", cf. Remark 2.16 
of |6|). Thus no such region Q." can exist. 

Proving both Case 4 and Case 5, therefore, involves showing that either O' or Q" must 
exist, unless w = in Q.. This can be done easily by reprising our arguments for ( ISll l- 
iS2\ . provided we know X. satisfies a maximum principle. Unfortunately, this will not in 
general hold true. However, if we let h(q, p) := h{q^, p), then by Theorem[T]we know that 
the maximum principle will hold precisely when Ai(£.) > 0. Lemma [T] then tells us that 
this will occur provided (12.20b and ( 12.21b hold. Making this assumption, the rest of the 
proof of the theorem is immediate. All that remains to show, therefore, is that iS3i implies 

dUQli-dliB- 

Let 6 1, 62 and b be defined as in Lemma [1] taking h(q,p) :- h{q'*^°,p) for Case 1, or 
h{q'^° + L,p) for Case 2. Using the equation satisfied by Ji, we can estimate 

sup \hpp\ < sup hhphqhqp - hphqq + gill - d(h))pphl -p{-p)hl\ 



ipp\ ^ aup ^z.,ip,iq,iqp — iipiiqq t — uyiijjfjpiip 

R R 

< 2Afal^ + Maq' + (sup ^1 + gT/max sup \pp\) 



where we take M := l|/!?llci(^)' which, in view of ( 12.7b . is precisely the same as in ( 11.131 ). 
Then, from the above estimate find, 

b, 62 < 4MAl + 2M^Al + 3Ao (sup + gi]^^ sup \pp\) 

+ 2M^ao'A^ + M^a„ 'A^ + Ma^^Al (sup |yS| + ^max sup \pp\) . 

The right-hand side, of course, we recognize as €2- Taking ei :- BCTmax sup |pp|, (12. 20b - 
(12. 21b follow immediately from (|S3b . This completes the theorem. □ 
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Remark. From the above proof we can see quite clearly why providing criteria for the 
symmetry of stratified waves is dramatically more difficult than for the homogeneous case: 
it is precisely the appearance of density variation that breaks the maximum principle struc- 
ture. Naturally, this leads us to wonder whether or not it is reasonable to expect that highly 
stratified waves are necessarily symmetric. While we are not prepared to ofl'er any conjec- 
tures on this question, it seems clear that, at the very least, making significant progress will 
require a new point of view. 
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